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Abstract: The resonant tunneling phenomenon is well understood in quantum mechanics. 
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1. Introduction 

In quantum mechanics (QM) the tunneling probability (or the transmission coefficient) of 
a particle incident on a barrier is typically exponentially suppressed. Somewhat surpris- 
ingly the addition of a second barrier can increase the tunneling probability for specific 
values of the particle's energy. This enhancement in the tunneling probability, known as 
resonant tunneling, is due to constructive interference between different quantum paths of 
the particle through the barriers. Under the right conditions, the tunneling probability can 
reach unity. This is a very well understood phenomenon in quantum mechanics [1]. The 
first experimental verification of this phenomenon was the observation of negative differen- 
tial resistance due to resonant tunneling in semiconductor heterostructures by [2]. In fact, 
this phenomenon has at least one industrial application in the form of resonant tunneling 
diodes [3]. 

Tunneling under a single barrier in quantum field theory (QFT) with a single scalar 
field is well understood, following the work of Coleman and others [4,5]. Despite some 
arguments given in [6,7], the issue of resonant tunneling in quantum field theory remains 
open [8,9]. Recently, Sarangi, Shiu and Shlaer suggested that the functional Schrodinger 
method should allow one to study this resonant tunneling phenomenon [10]. In this paper, 
we apply this approach to study resonant tunneling in QFT with a single scalar field. We 
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Figure 1: A one-dimensional potential in quantum mechanics with three local minima separated 
by two barriers. Consider an incoming particle from the left. The resonant tunneling effect can 
take place in the tunneling from A ^ C via B. With the appropriate energy for the particle, the 
tunneling probability or transmission coefficient may be as large as unity. 



show that resonant tunnehng in quantum field theory does occur and describe its properties. 
Following Coleman, we shall work in the thin-wall approximation. 

Before going into any details, it is useful to give an intuitive argument why some effect 
like resonant tunneling should happen in QFT. Consider the following tunneling process 
for a potential shown in Figure ||. Let the tunneling rate from A ^ B he Ta—tB = De~^ , 
while its tunneling probability Ta^b = Ke~'^ , which are taken to be exponentially small. 
The prefactor D or K is of order unity with the proper dimension. Here we shall focus on 
the exponential factor. Suppose the tunneling rate from B ^ C is given by Tb^Ci which is 
also exponentially suppressed. Both Ta^b and Tb^c are evaluated using standard WKB 
method. A naive WKB analysis will suggest that the tunneling from A ^ C is doubly 
exponentially suppressed, i.e., Ta^c ~ Ta^bTb^c- However, this is not correct. Consider 
the typical time, namely Ia^c^ it takes to go from A ^ C . It should be the sum of the 
time it takes to go from A ^ B plus the time it takes to go from B ^ C. Since the typical 
time is simply the tunneling (or decay) time, which is the inverse of the rate of tunneling, 
we have 

-^=tA^C = tA-.B+tB^C = Y^ + T^ (1-1) 
-L A^C i A^B i B-*C 



So it follows that 



Ta^c « J^-f"-^ (1.2) 

J^A^B + J^B^C 



which is clearly not doubly suppressed. For the special case where T = Ta~>b = Tb^c-, 
we see that Ta^c ~ T/2. 
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In QM, the exponential enhancement from Ta^c ~ T"^ to T/2 is due to the resonant 
tunnehng effect. At the resonances, 

Ta-.c ~ 1 while it is off resonances. For 

a generic incoming wavefunction, with a spread in energy eigenvalues covering one or 
more resonances, this resonant effect yields Ta^c ~ T/2, so the resonances typically 
dominate the tunneling process A ^ C. More generally, the relations ( |1.1D and ( |1.2D are 
reproduced [6,7]. Since the argument for ( |1.1D is general, it should apply in QFT as well 
as in QM. This suggests that some phenomenon like resonant tunneling must take place in 
QFT. The challenge is to find and understand it. 

The functional Schrodinger method was developed by Gervais and Sakita [11] and 
Bitar and Chang [12]. It starts with the idea that tunneling is dominated by the most 
probable escape path (MPEP) developed by Banks, Bender and Wu [13]. In QFT with a 
single scalar field 0, this path in the field space is described by 00 (x, A), where x stands 
for the spatial coordinates and A is a parameter that parametrizes the field configurations 
in the MPEP. In Coleman's Euclidean instanton approach, A is chosen to be the Euclidean 
time r, and the 0(4) symmetry of the instanton simplifies the analysis. On the other hand, 
the functional Schrodinger method allows one to make a different choice of A. In the lead- 
ing order WKB approximation, a generic choice leads one to a simple time-independent 
Schrodinger equation. It is not surprising that the resulting WKB formula for a single 
barrier tunneling process reproduces that of the Euclidean instanton approach. We shall 
present a step by step comparison so the equivalence of the two approaches is transparent. 
It is also obviously clear that the functional Schrodinger method is cumbersome by com- 
parison. However, this method has the great advantage of being immediately generalizable 
to the double (actually multiple) barrier case. The underlying reason is that the same real 
parameter A parametrizes both the under barrier (Euclidean time r) and the classically 
allowed (Minkowski time t) regions. 

For tunneling from vacuum A to vacuum C via the intermediate vacuum B in scalar 
QFT, we consider the simultaneous nucleation of two bubbles, where the outside bubble 
separates A from B and the inside one separates B from C. The functional Schrodinger 
method reduces this problem, in the leading WKB approximation, to a one-dimensional 
time-independent QM problem with A as the coordinate. The resulting double barrier 
potential in A, namely U{\), allows us to borrow the QM analysis to show the existence 
of resonant tunneling. In the case when both bubbles grow classically after nucleation, 
the tunneling process from AtoC will be completed. The tunneling rate is exponentially 
enhanced compared to the naive case. In the case where the inside bubble classically 
collapses back after its nucleation (because the bubble is too small for the difference in 
the vacuum energies to overcome the surface term due to the domain wall tension), only 
the tunneling from ^ to -B is completed. Still the tunneling rate from A to B can be 
exponentially enhanced. To draw a distinction between these two different phenomena, we 
call the second process catalyzed tunneling. For catalyzed tunneling, vacuum C plays the 
role of a catalyst. 

The rest of this paper is organized as follows. Resonant tunneling in quantum me- 
chanics is reviewed in Section ^. This review follows that in [1,6] . As we shall see, the 
functional Schrodinger method reduces the QFT problem to a QM problem, so the reso- 
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nant tunneling formalism in QM presented here goes over directly. In Section ^ we briefly 
review Coleman's Euclidean action approach, following [4]. In Section]^, we present the 
functional Schrodinger method. Here, the discussion follows closely that given by Bitar 
and Chang [12, 14] and we need only the leading order WKB approximation. For the 
single-barrier tunneling process, we see how Coleman's result is reproduced. In Section |5|, 
we discuss the double-barrier case. This is the main section of the paper. Here we find 
that resonant tunneling can occur in two different ways. It can enhance the tunneling from 
vacuum A to vacuum C via the intermediate vacuum B, or it can enhance the tunneling 
from ^ to S in the presence of C. Section ^ contains some remarks. The appendix contains 
a discussion on the particular ansatz for (p used in the main text. 



2. Review of Resonant Tunneling in Quantum Mechanics 

We first briefly review resonant tunneling in quantum mechanics. We consider a particle 
moving under the influence of a one-dimensional potential V{x) with three vacua shown in 
Figure Using the WKB approximation to solve the Schrodinger equation Hi/^ = Eij) for 
the wavefunction of the particle ip{x) gives the linearly independent solutions 

'iPl,r{x) ~ exp ^ ± i j dxk{x)^ (2.1) 



in the classically allowed region, where k{x) = ^/^{E — V{x)), and 



V'±(x) « ^=exp ( ± / dxK{x) ] (2.2) 

y/K{x) 



in the classically forbidden region, where k(x) = y ^{V{x) — E). A complete solution 
is given by il){x) = a^ipLix) + ORipn^x) in the classically allowed region and Tp{x) = 
a+V+(2;) + a-'ip-{x) in the classically forbidden region. To find the tunneling probability 
from A to B we need to determine the relationship between the coefficients ol,/? of the 
components ipL,R in vacuum A and the coefficients Pl,r in vacuum B. The WKB connection 
formulae give 




2l-i(e-e-i) e + e-i 




(2.3) 



where Q is given by 



~ 2 exp Q j ' dxy/2m{V{x) - E)^ , (2.4) 

and xi and X2 are the classical turning points. Setting (3l = 0, the tunneling probablity is 
given by 
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Since @ is typically exponentially large, Ta^b is exponentially small. 

The same analysis gives the tunneling probability from A to C, via B, as [1,6], 



where 



r^^c = 4j^(^GcI>+— J cos'w+i^- + -j sin^ Wj , (2.6) 
^> ~ 2exp Q y dxv^2m (y(x) - S)j (2.7) 



and 

1 /"^^ 



W = - dxJ2m(E -V(x)) , (2.8) 

with X3 and X4 the turning points on the barrier between B and C. 
If B has zero width, W = so Ta^c is very small, 

Ta^c - 46-2^-2 = Ta^bTb^c/^ (2.9) 

However, if W satisfies the quantization condition for the n^th bound states in B, 

W = {nB + l/2)-K (2.10) 

then cos = 0, and the tunneling probability approaches a small but not necessarily 
exponentially small value 

Ta^c = 9 (2.11) 

This is the resonance effect. If Ta^b and Tb^c are very different, we see that Ta-^c is 
given by the smaller of the ratios between Ta^b and Tb^c- Suppose Ta^b Tb^c- 
Following ( pl|) , we see that Ta^c 1 that is, the tunneling probability approaches unity. 



Notice that the existence of resonant tunneling effect here is independent of the detailed 
values of G, <I>, and W . 

The above phenomenon is easy to understand in the Feynman path integral formalism. 
A typical tunneling path starts at A and tunnels to -B. It bounces back and forth k times, 
where k = 0,1,2, ...00, before tunneling to C. When the Bohr-Sommerfeld quantization 



condition ( 2.10 ) is satisfied, all these paths interfere coherently, leading to the resulting 
resonant tunneling. On the other hand, if we raise the energy of the local minimum at B 
above the incoming energy E, then = and Ta^c ~ Ta^bTb^c which is typically 
doubly exponentially suppressed. 



3. Coleman's Euclidean Instanton Method 

Single-barrier tunneling in quantum field theory was studied by [4] using the Euclidean 
instanton method. For concreteness, we will focus on the (3 + l)-dimensional scalar field 
theory in Minkowski space described by the Lagrangian 

C = \^'-\{V<t^f-V{ct>) (3.1) 
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with an asymmetric double well potential 

Vicl)) = ^gi<P^-c'f-B{cP + c) (3.2) 

with B a small symmetry-breaking parameter. The potential at the false vacuum at (p = —c 
is zero, while the potential at the true vacuum at = +c is — e ~ —2Bc. The energy 
density difference between the two minima, namely e, is assumed to be small so that we 
may confine our analysis of this potential to the thin-wall regime. In the under-barrier 
(i.e., classically forbidden) region, one starts with the Euclidean action ^^(^^(t, x)) where 
T is the Euclidean time. Solving the resulting Euclidean equation of motion for cf) (with 
appropriate boundary conditions) and substituting it back into S'£'(0(r, x)) yields Se- In 
the semi-classical limit [4] showed that the tunneling rate per unit volume is 

r/V = Aex-p{-SE/h) , (3.3) 

where the sub exponential prefactor A studied in [5] will be unimportant for our purposes. 
The solution (/)(r, x, R) to the Euclidean equation of motion is the familiar 0(4)-symmetric 
domain-wall solution 

(I)dw{t, x, R) = -ctanh (^|(r - R)^ , (3.4) 

where /j, measures the (inverse) thickness of the domain wall, 

^i = (3.5) 

so a relatively large (i.e., thin wall) is assumed. Here 

r2 = |x|2 + r^ (3.6) 

and R is the radius of the bubble. The bubble wall sits at r = R, so cp = +c for r <C -R 
and = — c for r » R. That is, the bubble is surrounded by the false vacuum. It is useful 
to introduce the tension of the domain wall, 

S, = £ #v/2F(^ ^ £ rf<^yP^r^ = l^c (3.7) 

so the Euclidean action of this solution is now given by 

Se = -^TT^-R^e + "^TT^R^Si . (3.8) 

where the first term is the four-volume times the energy density difference e while the 
second term is the contribution of the domain wall. Setting the variation of this Se to zero 
yields 

£ = -^irR^e + AttR^Si = . (3.9) 

So the action is stationary for 

i? = Ac = 351/6, (3.10) 
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which gives us 

Se = Y'^i'^c = (^•^^) 

What happens to the bubble after nucleation? The bubble will behave in a way to decrease 
the energy £, i.e., d£/dR < 0. It is easy to see that the bubble prefers to grow (classically) 
as long as 

R > 2XJ3 (3.12) 

which is the case here. So, once the bubble is created with radius Ac, the domain wall starts 
at rest and moves (classically) outwards, eventually attaining relativistic speed. Notice that 
the condition ( |3.9D is simply the classical energy conservation equation: at the moment right 
after bubble nucleation, the total energy £ of the bubble and its interior equals to that of 
the original false vacuum in the region, which is zero. 

4. Functional Schrodinger Method 

Now let us introduce the functional Schrodinger method and apply it to the single barrier 
tunneling process discussed above in the same scalar QFT. In the semi-classical regime, 
a discrete set of classical paths, namely the most probable escape paths (MPEP) in con- 
figuration space give the dominant contributions to the vacuum tunneling rate [11-13]. 
Essentially this approximation allows us to reduce an infinite-dimensional quantum field 
theory calculation to a one-dimensional quantum mechanical computation. The effects of 
nearby paths can be included systematically in an h expansion, and were calculated to 
0{h?) in [11,12], but will not be relevant for the rest of our analysis. 

The Hamiltonian for a scalar field </>(t,x) where x denotes the three spatial directions 

is 



H = l d'^(^^ + ^{V<t>r + V{<P)^ . (4.1) 



where V{(f)) is that given in ( |3.2| ). To quantize the field theory we use [(/)(x), 0(x')] = 
iM'^(x — x') to replace (j) with —ihS/6(j). This replacement allows us to write the time- 
independent functional Schrodinger equation as 

iJ^'((/>(x)) = S^'(0(x)) (4.2) 

where 



and the eigenvalue E is the energy of the system. As usual \I'((/)(x)) is the amplitude that 
gives a measure of the likelihood of the occurance of the field configuration 4'{x). 

With the ansatz ^{(j)) = Aexp(— ■|S'(0)) where A is constant, the functional Schrodinger 
equation (|4.2|) becomes 



d^x I — 



2 



h~6^ ¥ \ ) 



+ ^(v</,)2 + V{^)\ e^^W = i?e^^W . (4.4) 
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Expanding £'((/>) in powers of h, S'(</>) = 5(o)(</') + ^'S'(i)(</') + and comparing terms with 
equal powers of h the functional Schrodinger equation ( [4.4D yields 



d^x 



1 fSS^o){^) 



--E, 

= 0, 
etc. 



(4.5) 



The infinite set of nonlinear equations (4^) on an infinite-dimensional configuration space 
can be reduced to a one-dimensional equation in the leading approximation. For our 
purpose here, we shall focus on 5(o) and ignore the higher-order corrections 5(2), etc. 
The essential idea is that there is a trajectory in the configuration space of (/'(x), known 
as the most probable escape path (MPEP), perpendicular to which the variation of 5(o) 
vanishes, and along which the variation of 5(o) is nonvanishing. We use A to parametrize 
this path in the configuration space of (j){x), so the MPEP is A). This MPEP satisfies 

^S(o) _ ^f.. 0(1)0 



6S, 



(0) 



(x,A) 







Along the MPEP, we have 



dS, 



(0) 



dX 



3 a(/>o(x,A) 55'(o) 



(4.6) 



(4.7) 



so C(A) is determined and we have 



(0) 



95, 



l0o(x,A)= 



(0) 



dX 



d^x 



d4>o{x,X) 
dX 



(x,A) 



dX 



We now define the effective potential U{X) = U{(j)(x, A)), so 



u{H^,x)) 



d^xl i(V<^(x,A))2+T/(</>(x,A)) 



(4.8) 



(4.9) 



then the classically allowed regions have C/((/)o(x, A)) < E and the classically forbidden 
regions have C/((/)o(x, A)) > E. 

Using the zeroth-order equation in ( |4.5| ) with ( |4.8| ) we find the WKB equation 

\ -1 



d^x 



dX 



dX 



C/(0(x,A))-i?. 



(4.10) 



To find the WKB wavefunctional it is sometimes useful to rewrite ( [4.10 ) in terms of the 
path length ds in the configuration space. This path length is defined by 

2 



{dsY 



d^x 



g0(x,A) 
dX 



{dXYm{(P{x,X)) . 



(4.11) 
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That is, choosing s to parametrize the MPEP, ( 4.10| ) simphfies to 

1 rds^o)V 



ds 



U{(t>{^,s))-E . (4.12) 



Now let us consider first the classically forbidden region, U{(j){x,s)) > E. The solution to 

(EH 



IS 



5(0) = i dsV2[f/(0(x, s')) -E]=i jj' d\ V2[f/(0(x, A)) - E] . (4.13) 

where Xti and Xt2 are the turning points. Treating both ^((/)(x. A)) ( [4.1l|) and U{(p{x,X)) 
( |4.9| ) as functionals of (f), the Euler-Lagrange equation for i?i>(x. A) follows from setting the 
variation of 5(0) to zero. It turns out that the resulting equation of motion derived from 



(4.13) simplifies considerably if we choose r as the parameter where 



- = y/2[U{<t>{x,T))-E] . (4.14) 

With T as the parameter, setting the variation of ( [4. 13] ) equal to zero now yields 

a20(x,T) , ^2^, , 5F(0(x,r)) 



+ V^</)(x, r) = (4.15) 



Here r simply plays the role of Euclidean time, and (4.15) is simply the Euclidean equation 



of motion for (/)(x, r). Once we obtain the solution for (/)(x, r), we insert this solution into 
(4.13) to obtain the value Sq for 5(o). So we see that the functional Schrodinger method 



leads both to a determination of Sq along the MPEP and an equation that determines 
MPEP, namely, (/)o(x,r) itself. The subscript "0" indicates that it is the MPEP in the 
leading WKB approximation. The equation ( 4. 151) has 0(4) symmetry, so it reproduces 



the familiar 0(4) symmetric domain- wall solution ( |3.4[) in the thin- wall approximation (r 
is the four-dimensional radial coordinate, r"^ = t"^ + |xp), 

,^o(x, r) = -ctanh (^{r - A,)^ (4.16) 

after imposing the boundary conditions, (/>^— casr— >oo and (/»^-|-casr— >0 and 
^0(0' 
dr 



^^^^^ — 0. Recall that Ac is the critical radius of the bubble given by (3.10). 



For our purpose, it is now convenient to introduce the parameter A 

A = v'Af^ (4.17) 

which is the spatial radius of the bubble as shown in Figure Here 0o(x, A) = for 
A < 0. For Ac > A > and near the domain wall at r = y^^|x|^~-|-t^, we have r — Ac ~ 
±(|xp — A^)/(2Ac) ~ ±(|x| lb A)A/Ac. The corrections introduced by this approximation are 
exponentially suppressed far from the domain wall, so we can express the 0(4) symmetric 
solution (p^ as, 

0o(x, A) « -ctanh (^{\x\ " A)^) , (4.18) 
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Figure 2: Spatial slices of the familiar domain- wall solution ( 4.18| ) can be parameterized by the 
spatial radius of the bubble of true vacuum, A. Here A is the length of the horizontal dashed line. 
The radius Ac of the (bottom half) bubble is related to the Euclidean time r via A^ = + A^. The 
choice of the A parameter allows us to go smoothly from the classically forbidden region (below the 
X-axis) to the classically allowed region (above the x-axis). 



This MPEP solution is plotted in Figure ^. (Slight care is needed for A ~ 0, in which case, 
we simply go back to ( [4.161) .) For A < the solution is the false vacuum (/>o(x, A) = — c . 
As A > increases, quantum fluctuation tends to fluctuate towards the true vacuum. At 
A = Ac, the bubble is created, which then evolves classically for A > Ac. 

We have effectively reduced the WKB wavefunctional in the classically forbidden region 
to a WKB wavefunction which can be written using (4.13) as 



^((/)(x,A)) =^e*^°/'^ 



Aexp(-i 



d\^/2m{\)[U{\) - E] 



where m(A) is obtained by substituting the MPEP i;^o(x, A) (|4.18|) into m(0(x. A)), 



m(A) 



)(x,A) 



Ac 



(4.19) 



(4.20) 



Here m(A) is the effective mass, which is manifestly positive, and the second equality is 
obtained in the thin- wall approximation. It is also straightforward to evaluate U{\) by 
substituting MPEP (/)o(x,A) (|4l8|) into C/((/)(x,A)) (U), 

27r5i 



U{\) 



Ac 



(4.21) 



This approximation for U{X) is plotted in Figure H. The classically forbidden region for 
zero-energy tunneling is < A < Ac- As mentioned above, i;Ao(x, A) ( 4.18| ) needs correction 
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(i>m 



+c 










\ \ 

A=0 


-c 





Figure 3: The MPEP is a path through field configuration space parameterized by A. Here, x is 
the spatial radius, = — c is the false vacuum and ((> = +c is the true vacuum. For the MPEP 
given by ( 4.1§| ): A < corresponds to the false vacuum, Ac > A > corresponds to the formation 
of the nucleation bubble during the tunneling process, A = Ac corresponds to the completion of the 
nucleation of the bubble of true vacuum, and A > Ac corresponds to the classical growth of the 
bubble. 



for A ~ 0. The more accurate U{X) (the blue dotted curve) is also shown in Figure |5|. 
However, this difference is not important for our analysis here. 

Given m(A) ( [4.20| ) and U{\) ( |4.21| ), the tunneling probl em in QFT has been reduced 



to a one-dimensional time-independent QM problem. We can now perform the integral in 
the exponent in ( [4.19| ) and obtain, for E = 0, 

-iSo = ^4 = (4.22) 
4 e'^ 



which reproduces ( 3.11| ). (The factor of two difference is because, here, we are evaluating 



the exponent of the tunneling amplitude instead of the tunneling rate.) This result is 
completely expected since the expression ( 3.11| ) is obtained by integrating the Langrangian 



density with the 0(4) symmetric solution (|4.16| ) or ( |3.4| ) in polar coordinates while here 
we are performing the same integral, first in the spatial coordinates x and then in the A 
(or equivalently in r) coordinate. For an 0(4) symmetric solution, the latter approach is 
unnecessarily cumbersome. 

So, before moving on, let us give a preview of the advantage of the functional Schrodinger 
method. This will also shed light on the underlying physics. In the classically allowed re- 
gions, similar arguments lead to the following 5(o), 



5(o)(<A(x,A)) = / dAV2m(</.(x,A))[ii;-[/((/<(x,A))] (4.23) 
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UiX) 




X 



Figure 5: The effective tunneling potential U{X) ioi g = 1, c — 1.5, and B — 0.1. The black 
solid curve shows the approximate form of U{X) given in ( 4.21 ). The blue dotted curve shows the 
numerical result. The difference is not important for our analysis. 



Similar to the previous case, the Euler-Lagrange equation of motion for cj) simplifies if we 
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choose the parameter t such that 



r7s 

- = x/2[i?-C/(,/)(x,t))] . (4.24) 



so setting the variation of 5(o) with respect to 0(x, t) to zero leads to 

a20(x,t) , , ay(0(x,t)) 



9^2 



V'</)(x, t) + ' = , (4.25) 



where t is simply the normal time, and this is simply the equation of motion for 4'{x, t) in 
Minkowski space, t > 0. Now, instead of t, let us choose A as the parameter, where 



which gives us (A = dX/dt), 



A = VA2+72 (4.26) 



^ ^ (4.27) 



Vl-A2 



which is simply the Lorentz factor. Substituting this into the path ( 4.1§| ), we obtain, for 
A > Ac, 

0o(x,A) = -ctanh(^|(|x|-A)A^ = _ctanh , (4.28) 

so the classical path cpoi'^, ^) now describes an expanding nucleation bubble, as its phys- 
ical radius A increases towards the limiting speed A = 1, with the proper Lorentz factor 



automatically included, and the effective tension now given by S'l/v 1 — A^. 

So we see that (/>o(x. A) for real A spanning oo > A > — cxd works equally well for 
classically allowed as well as classically forbidden regions. In summary, (/)o(x. A) = — c for 
A < 0, when (p stays in the false vacuum. For Ac > A > 0, 0o(x, A) describes the "averaged" 
quantum fluctuation that corresponds to the MPEP for the tunneling process under the 
potential barrier. The |x| < A region has fluctuated to the true vacuum, which is separated 
from the false vacuum region by a domain wall at |x| = A. Finally, for A > Ac, (f)Q{x,X) 
describes the classical propagation of the nucleation bubble. The advantage here is that a 
single real parameter describes the whole system. The problem has been reduced to that 
of a time- independent one-dimensional (i.e., the A coordinate here) quantum mechanical 
system of a particle with position-dependent mass m(A) ( 4.201 ) and a potential U{X) (^ 



which has a barrier (Ac > A > 0) that separates the two classically allowed regions. In 
general, the position-dependent mass complicates the quantization of the position variable. 
However, at leading order in the WKB approximation, such a complication does not arise. 



5. Resonant Tunneling in Scalar Quantum Field Theory 

The above discussion introduces the functional Schrodinger method and its applications to 
the tunneling process at leading order in h. It reduces the tunneling process in an infinite- 
dimensional field configuration space to a one-dimensional quantum mechanical tunneling 
problem. The essential difference between tunneling in field theory discussed in Section ^ 
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and in quantum mechanics discussed in Section |2| is that in field theory one should first 
find the MPEP, namely i;^>o(x, A), and then obtain the effective potential U{X) ( [4 .91 ) and the 
effective mass m(A) ( f4.20| ). We have extended the MPEP to include regions where classical 
motion is allowed. 

5.1 Setup 

To examine resonant tunneling in QFT, let us consider the following potential shown in 
Figure 

T/u^_ /i5i(('/' + ci)'-cf)2-i?i,/.-2i?iCi cjxO 

\i<72((<A-C2)2-ci)2-i?2<^-2i?iCi ^>0 ^'-'^ 

where as before Bi and B2 are small. For this potential the false vacuum (vacuum A) 
at (p ^ —2ci has zero energy density, the intermediate vacuum (vacuum B) at = has 
an energy density — ei = —2BiCi and the true vacuum (vacuum C) at ~ 2c2 has an 
energy density — ei — £2 = —2BiCi — 2B2C2. We take both ei and £2 to be small so that 
the thin-wall approximation is valid. Similar to the single barrier case, we introduce the 

^(i) - 2, 



inverse thickness fij = y2gjcj and the tension Sf = ^fJ'jCj for each of the two domain 
walls: j = 1 for the outside bubble and j = 2 for the inside bubble. Here ri > r2. In the 
thin- wall approximation, (j) = —2ci for r ^ ri, (j) = for ri ^ r ^ r2, and (j) = +2c2 for 
r <C ^2. The six parameters of the potential gi^2, ci,2 and Bi^2 now become /Ui_2, ' and 
61^2- In the thin- wall approximation, the thicknesses l//ii,2 of the domain walls drop out, 
simplifying the discussion. 

We also assume that the 0(4)-symmetric solution provides the dominant contribution 
to the vacuum decay rate. We note that the inside (half-)bubble does not have to be 
concentric with the outside (half-)bubble as long as the centers of the two bubbles lie on 
the same spatial slice. As we shall see, the analysis will go through without change as long 
as the two bubble walls are separated far enough, that is, much more than the combined 
thicknesses 1/ fii + 1/ fi2- We expect the off-center bubble configurations to be subdominant 
if we include corrections to the thin- wall approximation. We focus here on the zero-energy 
(i.e., E = 0) case. 

5.2 Ansatz 

The MPEP involves (j) in the two under-the-barrier regions as well as the classically allowed 
region between them. In the under-the-barrier regions, we can solve for </>(x, r) using the 
Euclidean equation of motion ( [4.15 ), while in the classically allowed region, we can solve 



for (/>(x, i) using the equation of motion in Minkowski space (4.25). We can then convert 
them to 00 (x, A). 

However, in the thin-wall approximation, it is easier to simply write down the ansatz in 
the radial coordinate and then extract (/>(x. A) from it. Here the tunneling process involves 
two concentric bubbles: an outside bubble whose domain wall separates A (outside) from 
B and an inside bubble whose domain wall separates B from C (inside). The radii of the 
two bubbles are ri and r2 as shown in Figure 0. As a function of the four-dimensional 
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radial coordinate r, we have the MPEP (see Appendix A), 

(/)(r) = -ci tanh (^^{i" - n)^ - C2 tanh (^^{f - ''2)^ + C2 - ci (5.2) 

For appropriate ri and r2, this solves the Euclidean equation of motion ( 4. 151) (and the 
Lorentzian equation of motion ( 4.25| ) in the appropriate regions). However, here we shall 
use this (p ( ^.2]) as an ansatz to find the resonant tunneling condition. 

Now it is straightforward to extract i?!)(|x|,A) from ij^(r) given by ( |5.2[) , 



(/)o(|x|, A) = -ci tanh ^y^d^l - A)^ - C2 tanh ^Y^*-'^' ~ '^''*) + - ci (5.3) 

where we use the same reparametrization as in the single-barrier case. Here A is the value 
of A at which the inside bubble has zero spatial extent, 

A2 = - r| (5.4) 

and as long as both bubbles are expanding 

I otherwise. 

This is shown in Figure 0. The equation ( ^.3]) also implies that (/)o(|x|, A) = — 2ci for A < 0. 
Note also that the sum of the second and third terms in ( ^.31 ) vanishes for A < A. Sub- 
stituting this MPEP i?i>o(|x|, A) given by ( ^.3D into ( [1.9| ) now yields, after a straightforward 
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T 




Figure 7: The tunneling process from A to C via B leads to the formation of two bubbles: the 
outside bubble separates B from A and the inside bubble separates C from B. They are drawn as 
concentric bubbles here, though this is not the case in general. Here A is the length of the horizontal 
dashed (red) line. We have ri > A > r2. 



calculation, the effective tunneling potential U{X) = U{<j)o{\x\, X)), 

UW = 2vr5« (A + ^) _ ^,,A3 + 2vr^f ) + ^) (V)^ - ^^.(A')^ (5.6) 

For appropriate parameter choices, we see that (^^) has four zeros as shown in Figure |8|. 
If the negative part of U{X) (i.e., the classically allowed region) is not too deep, Aa ^ A, 
and Xb ^ Aic = 3S[^^ /ei. We may take U{X) =0 for A < 0. The discontinuity in the 
derivative of U (A) at A will be smoothed when the thickness of the bubble wall is taken 
into account. 

It is also straightforward to evaluate the effective mass m(A) defined by ( |4.20| ) using 
(j)o{\x\,X) of dU), now given by 

^X' + ^{X')Ax. (5.7) 
ri r2X' ) 

Note that, as expected, m(A) > 0. Now we have a time-independent one-dimensional (with 
A as its coordinate) QM problem with the double-barrier potential C/(A) (|5.6|) and mass 
m(A) (^3), which is illustrated in Figure p. 

5.3 Constraints 

We see that the existence of the classically allowed region B' in \J (A) will lead to resonant 
tunneling. We would like to see what properties of potential T^((/') ( ^.1| ) will yield resonant 
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UiX) 




X 



Figure 8: A typical effective tunneling potential UiX) (5.6) for double tunneling. The classical 
turning points for zero-energy tunneling are 0, Ab, Aa, and ri. The discontinuity in the derivative 
of U{X) occurs at A = A. If the negative part of U{X) (i.e., the classically allowed region) is not too 



deep, As > Aic = 3S'P^/ei, and Aa > A. In appropriate units, the values plotted here are 



(1) 



^ = 1.4, ei = 0.25, £2 = 8.4 • 10"^ = 14 which gives n = 20 and A = 14.3 using (|5j) and 

(PD. 



tunneling. It is important to emphasize that the existence of the classically allowed region 
B' is not guaranteed. For E = the existence of a double-barrier U (A) potential requires 
four distinct classical turning points satisfying 

= U{0) = U{Xb) = U{Xa) = U{n) . (5.8) 

The radii of the two bubbles at the moment of nucleation are related via (| 



We may evaluate the Euclidean action Se for (p (5.2) in the thin-wall approximation. 
Minimization of Se with respect to ri and r2 separately will yield ri = Aic and r2 = X2c 
(see Appendix A). However, this bounce solution includes paths that passes through region 
B only once. Since resonant tunneling must include paths that bounce back and forth any 
number of times in the region B, this is not what we are seeking. Instead of finding the 
Se that includes multiple passes through B, we use the the functional Schrodinger method 
to reduce the problem to a one-dimensional time-independent QM problem, which is then 
readily solved for 5(o)- 

Once the simultaneous nucleation of the two bubbles is completed and just before they 
start to evolve classically, we are at A = ri, where U{ri) = 0, 



£(2) = U{n) = 47r(5f ) - \nei)rl + Ai,{sf^ - \r2e2)Tl = (5.9) 



This turns out to be the energy conservation condition as well. For the single-bubble case, 
the corresponding energy conservation condition ( |3.9| ) is equivalent to the minimization of 
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the action. For the double-bubble case, the total energy £(^2) of the two bubbles at the 
moment of creation (at A = ri) must vanish. If we treat the region between the bubbles 
classically during the nucleation process, then the energy of the inside bubble, that is, 
the second term in the above condition (|5.9| ) must vanish by itself (following from the 



condition (|3.9|) ), in which case the first term vanishes as well. That is, ri = Aic = 3S[^^ /ei 
(2) 

and r2 = X2C = 3S\ /e2- However, it is crucial that the classically allowed region receives a 
full quantum treatment. So we must treat the simultaneous nucleation of the two bubbles 
quantum mechanically and demand only if (2) = 0. The requirement that the Euclidean 
action be stationary in this case reproduces (5.9). See Appendix A for some details. 



Note that the existence of a classically allowed region B' implies that, U{X) < for 



A > A > Xb- Following (5.6), we obtain 



> A| = (5.10) 
2ri - Ale 



from which it follows that (A^ = '''i — '''2) 

n > r2 (5.11) 

ri > As > Ale = 35fVei 
r2 < X2c = 3SP/e2 . 

The existence of a second classically forbidden region requires lim^^^^- dU/dX\\ < since 
U{ri) = is automatically satisfied. Equivalently 

2(5f)+5f^) <riei+r2e2 . (5.12) 
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1-2 c 



2 A, 



A«=A 



No classically allowed region 




No classically allowed region 



Catalyzed 
Tunneling 

Z*****^^"^^ Dominant 
contribution 



dA I'^i 



Only one classically 
forbidden region 



ri 



Figure 10: The allowed parameter region for resonant tunneling. In appropriate units, the values 



plotted here are s[^'' = 1, S\^' — 5, ei — 0.1, £2 = 0.3. The energy constraint (5^) constrains ri 
and r2 to lie on the black solid curve. The region to the left of the blue dotted curve labeled As = A 
is excluded since in this region B' does not exist. The region below the blue dotted curve labeled 
dU/dX\ri = is excluded since in this region there is only one barrier in U{X). The red dots and 
green dot satisfy the Bohr-Sommerfeld quantization condition ( 2.1C| ) in addition to satisfying the 
consistency conditions (5.E), ( 5.10| ), (5.11), and ( 5.12| ). The green dot with r2 ~ 0.54A2c provides the 
dominant contribution to the tunneling probability. Since this point lies below the orange dashed 
line, catalyzed tunneling occurs. The tunneling probability is —logT^^^ w 2.6 • 10^, which is 
exponentially enhanced compared to the naive single-barrier tunneling probability — log T'^^™'^' = 
Se ~ 1.3 • 10^ 



:(2) 



These constraints are illustrated in Figure 10. We see that the simultaneous nucleation of 
two bubbles can now be parameterised by a single parameter, say r2- When permitted, 
the sizes of the bubbles, namely ri and r2, will be such that the resonance condition 
W = {n + l/2)7r is satisfied. 

The width of the classically allowed region AXb' = Aa — Ab in ^7(A) decreases mono- 
tonically as r2 increases. The classically allowed region is a point when AA^/ = at some 
maximum value r'2,max of r2 (when ( ^.lOD is saturated). When r2 > r2,max; there is no 
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C/(A) U{X) 




A 



Figure 11: Two examples where resonant tunneling is absent. In appropriate units, the values 
plotted here are S'j;^'' = 1, S^'^ = 5, ei — 0.1, £2 — 0.3. The left plot shows the effective tunneling 
potential U{X) for r2 — O.8A2C > 7'2,max- No classically allowed region exists for this potential. 
The right plot shows the effective tunneling potential for r2 = 0.4A2c < r2.min- There is no second 
barrier for this potential. The discontinuities in the derivative of U{X) will be smoothed out when 
the thickness of the bubble walls are taken into account. 



classically allowed region in U{X). Similarly the width of the second classically forbidden 
region AAbarrier = ''^1 " increases monotonically as r2 increases. At some minimum value 
?'2,min of r2 ( ^.121) is Saturated, and the second barrier in U{X) becomes a single point. The 
condition that AAbarrier > is equivalent to the condition lim^^^- dU/dX\x < 0. Figure 



11 shows a typical effective tunneling potential in each of these two cases. 

After the simultaneous nucleation of the two bubbles quantum mechanically, the out- 
side bubble will grow so the tunneling out of vacuum A will complete. Now there are two 
possibilities for the inside bubble, depending on whether it has the critical size ( 3.121 ) to 



grow or not (note that the binding energy of the two bubbles is expected to be negligible): 

(1) A2C > ^2 > 2A2c/3, in which case the inside bubble will grow as well. Hence 
the tunneling from vacuum A to vacuum C will complete, although the outside bubble is 
expected to grow faster than the inside bubble. This is the analogue of resonant tunneling 
in quantum mechanics, so we refer to this tunneling process from AtoC via B as resonant 
tunneling when W = {n + 1/2) vr. 

(2) < r2 < 2A2c/3, in which case the inside bubble will collapse after nucleation, 
while the outside bubble will grow. In this case, the tunneling from A to B will complete. 
At a later time, tunneling from B to C will take place via a normal tunneling process. In 
this process, the presence of vacuum C can increase the tunneling rate from A to B hy an 
exponential factor compared to the naive rate given by ( p. 11 ). We refer to this tunneling 



process from ^ to S in the presence of vacuum C as assisted or catalyzed tunneling, since 
C plays the role of a catalyst. Note that in this region (|5.5| ) is modified for A > ri: 

' VA2 - A2 A < A < ri 

y/r'i +r'i-X^ ri < A < ^/r'f + r| (5.13) 
otherwise. 



A' 



rl 



The inside bubble shrinks for A > ri and disappears entirely when A = \/rf + 
5.4 Tunneling Probability 

Explicitly, as shown in Figure]^, the classical turning points (with E = 0) are 0, Xb, X\ 
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and ri. Then the single-barrier tunneling probability is given by 



where now 

e = 2 exp 
In the thin-wall limit with E = 



d\yj2m{\)U{\) 



e 



3 r,(l) 



In ( -T 1 = -\%S, 



(5.14) 
(5.15) 

(5.16) 



As before, the tunneling probability (2^) calculated using the functional Schrodinger equa- 
tion, with B given by ( 5.16| ), agrees with the result of the Euclidean instanton method, 
exp(— 5£;//i), where Se is given by ( ^.11 ). 

The tunneling probability from vacuum A' to vacuum C via the intermediate vacuum 
B' , is given by (|2.6| ) where now W and $ are given by 

Aa 



W 



dXy/2m{X){-U{X)) 



Xb 
S^A 



Xj,-s[^hB log 



Aa + 



\2 



X2 -, 



Xb 



- XI 



^f^lelog 



X2 



(5.17) 



with the classical turning points shown in Figure ^. The third equality in ( 5.17] ) is valid if 
the classically allowed region is shallow. In this approximation, we also have, with Aa ~ A, 



In 



ri 



Aa 

Se 



dXyj2m{X)U{X) 



(5.18) 



vr 
T 



where Se ( |A.6D is given in the appendix. Since A, ri and r2 are related by (5^) and (p^), 
we may consider 0(r2), W{r2) and $(r2) as functions of r2 only. 

The bubble sizes are dominated by the ones that satisfy the resonance condition ( 2.1C| ) , 
i.e., W = [ub + l/2)7r for the n^th resonance. With r2 satisfying this condition and the 



constraints shown in Figure 10, the resulting tunneling probability is now given by ( 2.11 ): 

4 



(e/$ + «>/e)2 

which can approach unity for suitably chosen potential (^.1|). 

Next let us consider catalyzed tunneling. This is the case when the inside bubble 
classically re-collapses after its creation. The normal probability T^^b has a bounce value 



Se smaller than that given by (|5.16|) , or 



> 4/6^ 



(5.19) 
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where Q is given by (111). The presence of vacuum C can lead to an enhanced tunneUng 
probabiUty, 

4 

rpres ^ 

^■"^ ^ (e/$ + $/e)2 

which can be substantially bigger than T^^^^^ if ~ <I>. On the other hand, the catalytic 
effect is negligible if ^ is exponentially too big or too small when compared to 0. 

5.5 Generic Situation 

For large and <I>, so that the penetration through the barriers is strongly suppressed, 
the tunneling probability has sharp narrow resonance peaks at the values in ( 2.10D . If 



we allow the possibility of a non-zero inital energy E that is small compared to all other 
relevant mass scales, we introduce an extra variable without introducing any additional 
constraints, although (5.8) and ( |5.9| ) will be slightly modified. Treating the resonance 



shape as a function of energy E, the resonance has a width T^- Expanding around the 
resonance at E = E^, we have 



w ■ ■ 

cos W 



dE 
and 



(E-Er), sinTy = l (5.20) 



^E = —^{^ + ^) (5.22) 



so this yields, for large and ^, 

"~ 0$(|f) ' 
Next, let the separation between neighboring resonances be AE, where 

dE ) 

Then a good estimate of the probability of hitting a resonance is given by 

- ^' = If is (I + 1) = 2^ (^-^ + f^-^^' 

We see that the probability of hitting a resonance is given by the larger of the two decay 
probabilities, Ta^b or Tb^Ci and the average tunneling probability is given by 

< Ta^c >= P{A ^ C)Ta-.c ~ ^-^-^^^-^ (5.25) 

which is essentially given by the smaller of the two tunneling probabilities. This is a 
derivation of Eq.( |1.2D . Following the argument for (|1.1| ), the generalization to tunneling 
with multiple barriers is straightforward. 
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6. Remarks 



Our results do not conflict with the no-go theorem of [8] since the assumptions are inap- 
plicable as anticipated by [10]. (Note that in [8] the term MPEP is used exclusively to 
refer to the path in the classically forbidden region; our more inclusive definition will be 
used in the discussion here.) In particular one of the assumptions is that the MPEP is 
everywhere stationary = or = at the boundary between the classically allowed 
region and the classically forbidden region, i.e. for A = and A = Aa. This condition 
is clearly violated by our MPEP (|5.3| ) which is only stationary everywhere for A < and 
A = ri. The exact solution ( ^.2| ) also violates this condition. 

In [9] it was shown how a bubble of vacuum A surrounded by vacuum B could produce 
a bubble of vacuum C with probability of order unity under certain conditions. It was 
assumed that the vacuum energy density of vacuum C was greater than the vacuum energy 
density of vacuum A or vacuum B. The inhomogeneous initial state violates one of the 
conditions of the no-go theorem. The physics is quite different because the asymptotic 
false vacuum is intermediate in field space. Since our analysis applies only in the double 
thin-wall approximation, it is possible that the oscillons or its quantized version may play 
an important role in a different setup. 

We apply the functional Schrodinger method to show how resonant tunneling takes 
place in quantum field theory with a single scalar field. Our analysis is carried out in the 
double thin-wall approximation. The double-barrier potential problem in QFT is reduced 
to a double-barrier potential problem in a time-independent one-dimensional QM problem, 
so the quantum mechanical analysis can be applied. 

The relevance of resonant tunneling is obvious if the potential has many local minima, 
as is the case of the cosmic landscape in string theory. So resonant tunneling in the presence 
of gravity is a very important question to be addressed. 

What happens if the conditions (^^), ( ^.101 ), ( 5.11| ), and ( 5.12 ) cannot be satisfied? 



Generically, the double thin-wall approximation breaks down and a more careful analysis 



is needed. Based on the analysis of (1J-), we are led to believe that the resonant tunneling 
phenomenon will continue to happen. So it is interesting to study more general cases to 
obtain a complete picture. 

A. Euclidean Bounce Solution 

Recall that, for a 0(4)-invariant bounce, with r'^ = + |xp, the Euclidean equation of 
motion becomes 

£4> ^ 3 # ^ dVEjc^) 
dr"^ r dr d<j) 
where the Euclidean potential Ve{4>) = —^(0) and 

lim (t){r) = (I)A = -2ci (A.2) 

r— >oo 

Also, 

Ftrh 

(A.3) 



dr 



r=0 
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To demonstrate that a solution always exists, we simply follow Coleman's argument [4] 
that there are undershoot solutions as well as overshoot solutions. Continuity then implies 
the existence of the solution. If we start with (j){0) < (ps where (ps is close to but on the left 
of 4>c = 2c2, where Ve(0s) = VE{(t>B) = Ve{0) (that is Ve((/>(0)) < Ve((/)b)), then because 
of the damping term, we shall undershoot. For overshoot, we start at (pt very close to (pC: 
so we can linearize ( [A.l[ ), 

+ z:jz-f^ ]i^-<f>s) = o (A.4) 



^ dr'^ r dr 

where /x^ = V^{cpc)- The solution can be expressed in terms of a Bessel function, 

(P{r) -<Ps = 2[^s - (pc]h{f^r)/^r (A.5) 

where <ps — (pc < so is moving to the left. For (pg sufficiently close to </>c, we can arrange 
for (p to stay arbitrarily close to (pc for arbitrarily large r. For large enough r, the damping 
term is negligible. Without damping, </) can simply move beyond (pA- 

To satisfy the boundary conditions ( |A.2| , A.3 ), we expect a unique solution. We may 



evaluate the Euclidean action Se and then minimize it, which is equivalent to solving the 
above system ( [A.l - A.3). In the thin-wall approximation, the solution will take the form 



given by (|5.2|). This yields 



^2 , , ^2 



Se = 2^\lS^^ - \eM - 4) + 2^\lsf^ " y(^i + ^2)r| (A.6) 

A simple minimization of the resulting Se with respect to ri and r2 separately will yield 

5Se 



6rj 



0^rj = Xjc (A.7) 



However, this bounce solution includes paths that passes through region B only once. It is 
crucial to recognize that possible resonant tunneling must include, in terms of Feynmann's 
path integral formalism, paths that bounce back and forth any number of times in the 
region B. So, we are not interested in minimizing the action with a single pass through B. 
Instead, we should leave ri and r2 free at this moment. Once the simultaneous nucleation 
of the two bubbles is completed and just before they start to evolve classically, energy 
conservation demands the condition ( |5.9D . For the single-bubble case, the corresponding 
energy conservation condition (^]^) is equivalent to the minimization of the action. For 
the double-bubble case, we can obtain the condition (|5.9D by varying Se with respect to 
dri = dr2 ■ 

Instead of finding Se that includes multiple passes through B, we use the functional 
Schrodinger method to reduce the problem to a one-dimensional time-independent QM 
problem, which is then readily solved. Here, the energy conservation condition ( ^.S| ) is 
simply the turning point of U{X). 
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